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We use observational data from Type la Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), 
and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis 
(BBN), to constrain the running parameter A of Hofava-Lifshitz gravity, which determines the flow 
between the Ultra- Violet and the Infra-Red. We consider both the detailed and non-detailed balance 
versions of the gravitational sector, and we include the matter and radiation sectors. Allowing for 
variation of all the parameters of the theory, we construct the likelihood contours and we conclude 
that in la confidence A is restricted to |A — 1| < 0.02, while its best fit value is |A6./ — 1| « 0.002. 
Although this observational analysis restricts the running parameter A very close to its IR value 1, 
it does not enlighten the discussion about the theory's possible conceptual and theoretical problems. 
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I. INTRODUCTION 

Hofava recently proposed a power-counting renormal- 
izable. Ultra- Violet (UV) complete theory of gravity ilr 
Although presenting an Infra-Red (IR) fixed point, 
namely General Relativity, in the UV the theory pos- 
sesses a fixed point with an anisotropic, Lifshitz scaling 
between time and space. Since then there has been a 
significant progress in examining the properties of the 
theory itself, including various extensions of the original 
basic version [5l-l30| . Additionally, application of Hofava- 
Lifshitz gravity as a cosmological framework gives rise to 
Hofava-Lifshitz cosmology (sil . |32| , and in particular one 
can study specific solution subclasses [33- 47], the phase- 
space behavior [48l - [5^ , the gravitational wave production 
[S^ - Issj . the perturbation spectrum [59l - l68j . the matter 
bounce [69l-[73|. the black hole properties [75l - [89| . the 
dark energy pheno meno logy ,90- 95], the astrophysical 
phenome nology 96l - llOCll |, the thermodynamic properties 
[lOlMlOSj etc. However, despite this extended research, 
there are still many ambiguities if Hofava-Lifshitz gravity 
is reliable and capable of a successful description of the 
gravitational background of our world, as well as of th e 
cosm ological behavior of the universe 

[li- 

Although the discussion about the foundations and the 
possible conceptual and phenomenological problems of 
Hofava-Lifshitz gravity and cosmology is still open in 
the literature, it is worthy to examine the possible con- 
straints that observational and cosmological data could 
impose on the parameters of the scenario (of the basic 
as well as of the various extended versions). In such 
an investigation a reasonable assumption is to set the 
running parameter A to its IR value 1, since all observa- 



tions lie deep inside the IR. Thus, under this assumption 
in [96l - [98| the authors used Solar System observations 
in order to constrain some of the remaining model pa- 
rameters of the basic version of Hofav a-Lif shitz cosmol- 
ogy under detailed balance, while in (lOOj | the analysis 
was performed including a soft detailed-balance break- 
ing. Similarly, in [lo3 we used cosmological observations 
(Type la Supernovae (SNIa), Baryon Acoustic Oscilla- 
tions (BAO) and Cosmic Microwave Background (CMB) 
ones, together with Big Bang Nucleosynthesis conditions) 
in order to impose complete constraints on all the param- 
eters of the basic version of Hofava-Lifshitz cosmology 
and construct the corresponding contour plots, with or 
without the detailed-balance condition. 

Although setting A to its IR value is a first and reason- 
able assumption, one could go beyond it, and examine the 
observational constraints that the data could impose on 
A itself. However, allowing A varying, that is preserving 
the Lorentz invariance breaking (which is restored in the 
exact IR value A = 1), one must take into account that 
in theories with Lorentz invariance breaking the "gravi- 
tational" Newton's constant Ggrav, that is the one that 
is present in the gravitational action, does not coincide 
with the "cosmological" Newton's constant Gcosmo, that 
is the one that is present in Friedmann equations 110]. 
Thus, in the case of Hofava-Lifshitz gravity one could 
use the deviation between Ggrav and G^r 



as it is con- 
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strained by measurements of the primordial abundance of 
He'* [llOl |. in order to extract an up per b ound on |A — 1|. 
This approach was followed in [20, Il06l | with the result 
< |A— 1| < 0.1. However, it is obvious that such an ap- 
proach can only provide a crude upper bound on |A — 1|, 
since it considers that all the other parameters of the the- 
ory remain constant. The correct approach should be to 
perform a systematic investigation, allowing for simulta- 
neous variations of all model parameters, and constrain 
all of them using observations. 

In the present work we are interested in performing 
such an holistic observational constraining of all the pa- 
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rameters of the basic version of Hof ava-Lifshitz cosmol- 
ogy, and especially of the running parameter A, using 
SNIa, BAO and CMB cosmological observations in or- 
der to construct the corresponding probability contour- 
plots. Furthermore, in order to be general and model- 
independent, we perform our analysis with and without 
the detailed-balance condition. The plan of the work is 
the following: In section |TT] we present the basic ingredi- 
ents of Hofava-Lifshitz cosmology, extracting the Fried- 
mann equations, and describing the dark matter and dark 
energy dynamics. In section IIIII we constrain both the 
detailed-balance and the beyond-detailed-balance formu- 
lations using cosmological observations, and we present 
the corresponding likelihood contours. Finally, in section 
IIVI we summarize the obtained results. 



II. HORAVA-LIFSHITZ COSMOLOGY 

In this section we briefly review the scenario where 
the cosmological evolution is governed by Hofava-Lifshitz 
gravity [Sl, 32]. The dynamical variables are the lapse 
and shift functions, N and Ni respectively, and the spa- 
tial metric gij (roman letters indicate spatial indices) . In 
terms of these fields the full metric is written as: 

ds^ ^ --N'^df + g^Jidx' + N'dt){dx^ + N^dt), (1) 

where indices are raised and lowered using gij. The 
scaling transformation of the coordinates reads: t — > 
l^t and lx\ 



A. Detailed Balance 

The gravitational action is decomposed into a kinetic 
and a potential part as S'g = / dtd^Xy/gN{£K+^v)- The 
assumption of detailed balance 0] reduces the possible 
terms in the Lagrangian, and it allows for a quantum 
inheritance principle [ij, since the {D + l)-dimensional 
theory acquires the renormalization properties of the D- 
dimensional one. Under the detailed balance condition 
the full action of Hofava-Lifshitz gravity is given by 



So 



dtd^x^N I -^{K^jK'^ - XK^) 



where 



8(3A- 1) 



+ KR- 3A' 



is the extrinsic curvature and 



V9 



(2) 



(3) 



(4) 



the Cotton tensor, and the covariant derivatives are de- 
fined with respect to the spatial metric gij. e^^^ is the 
totally antisymmetric unit tensor, A is a dimensionless 
constant and the variables k, w and /i are constants with 
mass dimensions —1, and 1, respectively. Finally, we 
mention that in action ([2|) we have already performed the 
usual analytic continuation of the parameters and w of 
the original version of Hofava-Lifshitz gravity, since such 
a procedure is requ ired in order to obtain a realistic cos- 
mology [H m il!, [lOl| (although it could fatally affect 
the gravitational theory itself). Therefore, in the present 
work A is a positive constant, which as usual is related 
to the cosmological constant in the IR limit. 

In order to add the matter component (including both 
dark and baryonic matter) in the theory one can follow 
two equivalent approaches. The first is to introduce a 
scalar field [sil |32{ and thus attribute to dark matter a 
dynamical behavior, with its energy density p„i and pres- 
sure Pm defined through the field kinetic and potential en- 
ergy. Although such an approach is theoretically robust, 
it is not suitable from the phenomenological point of view 
since it requires specially-designed matter-potentials in 
order to acquire an almost constant matter equation-of- 
state parameter {wm = Pm/ pm) as it is suggested by 
observations. In the second approach one adds a cos- 
mological stress-energy tensor to the gravitational field 
equations, by demanding to recover the usual general 
relativity formulation in the low-energy limit [l^, 0, . 
Thus, this matter-tensor is a hydrodynamical approxima- 
tion with Pm andpm (or pm and Wm) as parameters. Sim- 
ilarly, one can additionally include the standard-model- 
radiation component (corresponding to photons and neu- 
trinos), with the additional parameters pr and Pr (or pr 
and Wr)- Such an approach, although not fundamental, 
is better for a phenomenological analysis, such the one 
performed in this work. 

In order to investigate cosmological frameworks, we 
impose the projectability condition JJJ and we use an 
FRW metric 



iV = 1 , g,j^a^{t)-f, 



with 



jijdx^dx-' 



dr^ 



l-Kr^ 



r^dfll , 



(5) 



(6) 



where K <,=,> corresponding to open, flat, and 
closed universe respectively (we have adopted the con- 
vention of taking the scale factor a{t) to be dimensionless 
and the curvature constant K to have mass dimension 2). 
By varying N and gij, we extract the Friedmann equa- 
tions: 



K 



6(3A- 1) 
6(3A- 1) 
8(3A- l)2a2 



(pm + Pr 



8(3A-l)a4 8(3A-1) 



(7) 
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^+2^ - "4(3A-1) 



Yy (wmPm + WrPr^ - 



4(3A- 1) 

16(3A- l)2a2 



8(3A-l)a4 8(3A-1) 



(8) 



where H = - is the Hubble parameter. As usual, pm 
(dark plus baryonic matter) follows the standard evolu- 
tion equation 

p,n + SH{p,n + Pm) = 0, (9) 

while pr (standard-model radiation) follows 

Pr+3H{pr+Pr) ^0. (10) 

Lastly, concerning the dark-energy sector we can define 



Pde 



8(3A-l)a4 8(3A-1) 



= 8(3A- l)a4 " 8(3A-1)' 



(11) 



(12) 



The term proportional to a^^ is the usual "dark radia- 
tion term" , present in Hofava-Lifshitz cosmology [3TI,[33|. 
while the constant term is just the explicit cosmological 
constant. Therefore, in expressions (fTT |) . (|12p we have 
defined the energy density and pressure for the effec- 
tive dark energy, which incorporates the aforementioned 
contributions. Finally, note that using (|11I) . ((T2|) it is 
straightforward to show that these dark energy quantities 
satisfy the standard evolution equation: 



Pde + iH{pDE + Pde) = 0. 



(13) 



Using the above definitions, we can re- write the Fried- 
mann equations (I?]),® in the standard form: 



6(3A- 1) 



Pr + Pde 



8(3A- l)2a2 



(14) 



Pm+Pr+PDE 



k'^P^KK 



4(3A - 1) r '" ' ' "^^J 16(3A - l)2a2 ' 

(15) 

Therefore, if we require these expressions to coincide with 
the standard Friedmann equations, in units where c — 1 
we set [si [si: 



Gcosn 

8(3A- 1) 



167r(3A - 1) 



= 1, 



(16) 



where Gcosmo is the "cosmological" Newton's constant. 
Note that as we said in the Introduction, in theories 



with Lorentz invariance breaking Gcosmo does not coin- 
cide with the "gravitational" Newt on's constant Gg^av, 
unless Lorentz invariance is restored |llCl| . For complete- 
ness we mention that in our case 



G, 



grav 



327r' 



(17) 



as it can be straightforwardly read from the action ([2]) 
(our defi nitions of Gcosmo, Ggrav coincide with those of 
|20l Il06| ). Thus, it becomes obvious that in the IR 
(A = 1), where Lorentz invariance is restored, Gcosmo 
and Ggrav coincide. 



B. Beyond Detailed Balance 

The above formulation of Hofava-Lifshitz cosmology 
has been performed under the imposition of the detailed- 
balance condition. However, in the literature there is a 
discussion whether this condition leads to reliable results 
or if it is able to reveal the full information of Hofava- 
Lifshitz gravity [HI, HI]- Thus, one should study also the 
Friedmann equations in the case where detailed balance 
is relaxed. In such a case one can in general write [TT[ . 



2ao 



+ 



(3A - 1) 
2 

(3A^ 

0-2 



6 
K 



6a6 



3(3A- 1) a2 



(18) 



H 



3cto 



(3A 
3 



(3A-1) 
a'2 



6 



WrPr 



18a4 



6a6 



K 



6(3A- 1) o2 



(19) 



where ctq = '*^/12, and the constants Ui are arbitrary 
(with 0-2 being negative). Note that one could absorb the 
factor of 6 in redefined parameters, but we prefer to keep 
it in order to coincide with the notation of [ij, [4^ . As 
we observe, the effect of the detailed-balance relaxation 
is the decoupling of the coefficients, together with the 
appearance of a term proportional to . In this case the 
corresponding quantities for dark energy are generalized 
to 



Pde\ 



Pde\ 



6 



11 

' 6 



6a4 
18a4 



6a6 

CFiK 



(20) 
(21) 



Again, it is easy to show that 



iH{p 



DE\ 



+ Pl5i^l„„„-aJ=0. (22) 
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Finally, if we force p8)) . (|19p to coincide with the standard 
Friedmann equations, we result to: 



G, 



6ao 



87r(3A-l) 
(72 = -3(3A- 1), 



(23) 



while in this case the "gravitational" Newton's constant 



Ggrav reads [13|: 



IGtt 



(24) 



III. OBSERVATIONAL CONSTRAINTS 

Having presented the cosmological equations of a uni- 
verse governed by Hoi^ava-Lifshitz gravity, both with and 
without the detailed-balance condition, we now proceed 
to study the observational constraints on the model pa- 
rameters. This is performed in the following two subsec- 
tions, for the detailed and non-detailed bala nce s cenarios 
separately. We mention that, contrary to |l09l |. in this 
work we allow the running parameter A to vary too, and 
in order to be general enough we do not use any the- 
oretical argument to restrict it in any specific interval, 
handling it as completely free. 



A. Constraints on Detailed-Balance scenario 

We work in the usual units suitable for observational 
comparisons, namely setting 87rGgrav — 1 (we have al- 
ready set c = 1 in order to obtain ((16])). This allows 
us to reduce the parameter space, since in this case P7)) 
gives 



= 4, 



(25) 



and thus (|T6)l lead to: 

^cosn 



47r(3A - 1) 
= i^^. 



(26) 



Inserting these relations into Friedmann equation ([7]) we 
obtain 



3(3A- 1) 







3A\ 









(27) 



In terms of the usual density parameters {flm = 
p,n/i3H^), = ~K/{H^a^), rir = Pr/i3H^)) this ex- 
pression becomes: 

^-(3A^(^™+"^)-"--i^(2l^ + t)- 



Applying this relation at present time and setting the 
current scale factor oq = 1 we obtain: 



1 - 



(3A-1) 



n 



rnO 



QrO ) ^ ^KO — tt2 



1 



2A 



(29) 

where a 0-subscript denotes the present value of the cor- 
responding quantity. Note that flmo includes contribu- 
tions from both baryons Vlbo as well as dark matter JIdmo- 
In order to proceed to the elaboration of observa- 
tional data, we consider as usual the matter (dark plus 
baryonic) component to be dust, that is Wm ~ 0, and 
similarly for the standard-model radiation we consider 
Wr — 1/3, where both assumptions are valid in the epochs 
in which observations focus. Therefore, the correspond- 
ing evolution equations (|9)). (|10p give pm = Pmo/a^ and 
Pr = Pro I respectively. Finally, it proves convenient to 
use the redshift z as the independent variable instead of 
the scale factor (1 -f z = ao/a = 1/a). Inserting these 
into Friedmann equation (1271) we obtain 



r!„o(l + 2)' + r!.o(l + ^)' 



2( 2 

H(3A- 1) 

+nKo{l + zf+\uJ+^-^{l + zf]}, (30) 



where we have also introduced the dimensionless param- 
eter 



A 

Thus, the constraint ([29]) can be rewritten as 



(3A 



— (flmO + ^rO^ + ^KO + 1^ + 



f72 



(31) 



(32) 



We remind that the term O^^q/ (4a;) is the coeflScient of 
the dark radiation term, which is a characteristic feature 
of the Hofava-Lifshitz gravitational background. Since 
this dark radiation component has been present also dur- 
ing the time of nucleosynthesis, it is subject to bounds 
from Big Bang Nucleosynthesis (BBN ). As discussed in 
more details in the Appendix of 109], if the upper limit 
on the total amount of dark radiation allowed during 
BBN is expressed throug h the pa rameter AA''^ of the 
effective neutrino species |llll4ll4 |. then we obtain the 
following constraint : 



^ ^ 0.135AAf^f]^o- 

4a; 



(33) 



In this work, in order to ensure consistency wit h BB N, 
we adopt an upper limit of AA'^^ < 2.0 following |ll3l |. 



In most studies of da rk energy models it is customary 
to ignore curvature (e.g. [ll6l4l23j |). especially concerning 
observational constraints. This practice is well motivated 
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K 


K7(87rGgrav) 


(l/Hl) A 


(87rGg,avi?0) M 


A 




>0 


4 


(0, 1.46) 


(1.37, cx)) 


(0.98, 1.01) 


(0, 0.32) 


<0 


4 


(0, 1.46) 


(1.8, oo) 


(0.97, 1.01) 


(0, 0.68) 



TABLE I: 1(7 limits on the parameter values for the detailed- 
balance scenario, for positive and negative curvature. The 
cosmological parameters fimo, f^bo, ^ro and Ho have been 
marginalized over. 



since most inflationary scenarios predict a high degree 
of spatial flatness, and furthermore the CMB data im- 
pose stringent constraints on spatial flatness in the con- 
text of constant-w models (for example a combination 
of WMAP-(-BAO-|-SNIa data 115] provides the tight si- 
multaneous constraints —0.0179 < VIko < 0.0081 and 
-0.12 <l + w < 0.14, both at 95% confidence). 

However, it is important to keep in mind th at du e to 
degeneracies in the CMB power spectrum (see [l24j and 
references therein) , the limits on curvature depend on as- 
sumptions regarding the underlying dark energy scenario. 
For instance, if instead of a constant w one assumes a 
linearly varying w (that \s w [a) = wq + (1 — a) Wa) , the 
error on ^k u is of th e ord er of a f ew percent, that is 
much larger fl25l - [T27| (see fl28l - [T30} for the constraints 
on curva ture for different parameterizations) . The au- 
thors of |130| showed that for some models of dark en- 
ergy the constraint on the curvature is at the level of 
5% around a flat universe, whereas for others the data 
are consistent with an open universe with VIkq ^ 0.2. 
According to |l26t . geometrical tests such as the com- 
bination of the Hubble parameter H{z) and the angu- 
lar diameter distance Da{z), using (future) data up to 
sufficiently high redshifts z ~ 2, might be able to disen- 
tangle curvature from dark energy evolution, t hough not 
in a model-independent way. Furthermore, in | l3ll Il32j | 
the authors highlighted the pitfalls arising from ignor- 
ing curvature in studies of dynamical dark energy, and 
recommended to treat ^Iko as a free parameter to be 
fitted along with the other model parameters. Lastly, 
note that in the present work the spatial curvature plays 
a very crucial role, since Hofava-Lifshitz cosmology co- 
incides completely with ACDM if one ignores curvature 
[31I l32j . Therefore, and following the discussion above, 
we choose to treat VIko as a free parameter. 



In summary, the scenario at hand involves seven model 
parameters namely the cosmological parameters flmo, 
flbo, i^fcO; ^ro and Hq and the model parameters A, w and 
AN^, subject to constraint equations ([5^ and ([55)) . We 
marginalize over the cosmological parameters flmo, i^bOi 
riro and Hq. Of the four remaining parameters, only two 
are independent. We choose A and AN^ as our free pa- 
rameters. Once these are chosen, and for a given choice 
of curvature, ^ko and uj are immediately fixed from the 
constraint equations. In particular, lo can be determined 



L1 



L06 



L02 



0.98 



0.94 



0.9, 



K>0 




0.1 



0.2 
AN 



0.3 



0.4 



1.06 



1.02 




0.98 



0.94 



0.16 



FIG. 1: (Color online) Contour plots of X vs ANi, for nega- 
tive (K < Q) eurvature, in the detailed-balance scenario, un- 
der SNIa, BAO and CMB observational data. The remaining 
parameters have been marginalized over (see text). The yel- 
low (light) region is excluded at the 2a level, and the orange 
(darker) region is excluded at the la level. The red (darkest) 
region is not excluded at either confidence level. The white 
diamonds mark the best-fit point, and the A = 1 line is drawn 
for convenience. 



by eliminating VIkvi from relations ([32)) and ([33)l: 



cj-2sgn(f7Ko) \/Q.l'ibANynroUJ + Q.libAN^VLrQ 



3A- 1 



-1 = 0. 



(34) 



This fractional-order equation in w can have one or two 
roots. In the latter case we use the larger one, since our 
goal is to set an upper limit on uj (and hence on A). Vtxo 
can then be found from u using ([33]). 

In Fig. [1] we display the likelihood-contours for the 
free parameters A vs AN^, for both positive and negative 
curvature. All other parameters have been marginalized 
over. The details and the techniques of th e fitting proce- 
dure can be found in the Appendix of |l09j | . Additionally, 
in Table H] we summarize the la limits on the parameter 
values for the detailed-balance scenario. These bounds 
are in agreement with the corresponding ones of our pre- 
vious work I109I. Furthermore, as we observe, the data 
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constrain A to roughly A = l!to o2 the la level, while 
its best fit value is very close to 1 {Xb.f = 0.006) for both 
positive and negative curvature. 



B. Constraints on Beyond-Detailed-Balance 
scenario 



In units where SttG^ 



1 relation (p4)) gives 



CTo - 1/3. 



(35) 



Using this value and following the procedure of the pre- 
vious subsection, the Friedmann equation (I18|) can be 
written as 



H, 



«{ 



(3A - 1) 



+r!A'o (1 + zf + 

^ (3A^ K + (1 + -"^^ +iOi{l + zf^y, (36) 

where we have introduced the dimensionless parameters 
wi, wa and uj^, related to the model parameters cti, 
and (74 through: 



UJl = 

CJ3 = 
074 = 



6 

g"4^JfO 

6 



(37) 



Additionally, we consider the combination W4 to be posi- 
tive, in order to ensure that the Hubble parameter is real 
for all redshifts. a;4 > is required also for the stability of 
the gravitational perturbations of the theory [13, [l^- For 
convenience we moreover assume (73 > 0, that is ws > 0. 

In summary, the present scenario involves the following 
parameters: the cosmological parameters Hq, flmo, flxoi 
^bo, f^T-O) and the model parameters A, wi, W3 and UJ4. 
Similarly to the detailed-balance section these are subject 
to two constraints. The first one arises from the Friedman 
equation at z = 0, which leads to 



(3A-1) 



Qko^I- (38) 



This constraint eliminates the parameter wi . The second 
constraint arises from BBN considerations. The term in- 
volving wa represents the usual dark-radiation compo- 
nent. In addition, the LiJ4-term represents a kination-like 
comp onent (a quintessence field dominated by kinetic en- 
ergy [l35l Il3(| V If ANt^ represents the BBN upper limit 
on the total energy density of the universe beyond stan- 
dard mod el co nstituents, then as we show in the Ap- 
pendix of [109*1 we acquire the follo wing constraint at the 
time of BBN (z = zbbn) [liMIl: 



It is clear that BBN imposes an extremely strong con- 
straint on CJ4, since its largest possible value (correspond- 
ing to W3 = 0) is ^ lO^^'*. Finally, wamax denotes the 
upper limit on cja- In the following, we use expression 
(p9| to eliminate UJ4. For convenience, instead of W3 we 
define a new parameter 



UJ3 



Wan 



(40) 



which has the interesting physical meaning of denoting 
the ratio of the energy density of the Hof ava-Lifshitz dark 
radiation to the total energy density of Hofava-Lifshitz 
dark radiation and kination-like components at the time 
of BBN. 

We use relation (p9| to eliminate ^4 in favor of a and 
AN^, and treat A, a, flRo and AN^ as our free param- 
eters, marginalizing over Hq, flmOj ^bo and flro- Using 
the combined SNIa-|-CMB+BAO data, we construct like- 
lihood contours for different combinations of the above 
parameters, which are presented in Fig. [21 In each case, 
the free parameters not included in the plot have been 
marginalized over. The details and the techniques of 
the c onstruction have been described in the Appendix of 
[l09j |. with the only differ ence being that in the present 
work WMAP 7- year data jl33l | and the more recent Con- 
stitution Supernovae dataset |134l | have been used. In 
Table [TT] we summarize the la limits on the parameter 
values for the beyond-detailed-balance scenario. 







a 


A 


(-0.01,0.01) 


(0,2) 


(0, 1) 


(0.98, 1.01) 



wa + a;4 (1 + zbbn) 



1^311 



0.135A7V^17ro- (39) 



TABLE II: 1(7 limits on the free parameters of the beyond- 
detailed-balance scenario. The cosmological parameters QmO, 
flbo, ^ro and Ho have been marginalized over. 



Furthermore, as we observe, in la confidence the run- 
ning parameter A of Hofava-Lifshitz gravity is restricted 
to the interval |A — 1| < 0.02, for the entire allowed range 
of oja (that is of era)- Finally, the best fit value for A 
restricts |A — 1| to much more smaller values, namely 
\\b.f - 1| w 0.002. 



IV. CONCLUSIONS 

In this work we have constrained the running parame- 
ter A of Hoi^ava-Lifshitz gravity using observational SNIa, 
BAO and CMB data as well as considerations from BBN. 
In order to be general enough we have not used any 
theoretical argument to a priori restrict A in any spe- 
cific interval, handling it as a completely free parameter. 
Additionally, we have performed our investigation under 
the detailed-balance condition, as well as under its relax- 
ation. Finally, we have included the matter and radiation 
sectors following the usual effective fluid approach. We 
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FIG. 2: (Color online) Contour plots of different pairs of free parameters in the beyond- detailed-balance scenario, under SNIa, 
BAO and CMB observational data. In each case the parameters not included in the plots have been marginalized over. Color 
scheme as in Fig. QJ 



stress that we have let all the parameters of the theory 
to vary, performing an overall and holistic investigation, 
and we have not followed the naive approach, that is to 
keep everything fixed and vary only A in order to match 
observations. 

As we showed, Hofava-Lifshitz cosmology, either with 
or without the detailed-balance condition, can be com- 
patible with observations. We constructed the likelihood- 
contours for the involved free parameters, and we found 
that in la level A is restricted to |A — 1| < 0.02. As ex- 
pected, these bounds are one order of magnitude tighter 
than the corresponding ones arising from the considera- 
tion of primordial He^-abundance me asure ments in order 
to extract a crude upper bound [13, Il06l | . Finally, con- 
cerning the best fit value for A we obtain |Af,./ — 1| « 0.006 
for the detailed balance and |Ab./ — 1| « 0.002 for the 
beyond-detailed-balance cases. The aforementioned fea- 
tures act as additional arguments in favor of the present 
holistic investigation, where all the parameters of the 
theory are allowed to vary, comparing to the naive or 
partial ones in which only A varies. Lastly, it is inter- 
esting to note that these results, arising from cosmo- 
logical observations, are relatively close to the prelim- 
inary parametrized post-Newtonian (PPN) estimations 
for Hofava-Lifshitz gravity (0 < |A — 1| < 4x 10~^ 



[137|), despite the fact that the Solar System measure- 
ments (that lie at the basis of PPN parameters 138l - ll4l| ) 
are much more accurate than the cosmological ones. 

In summary, as expected, we found that the value of 
the running parameter A of Hofava-Lifshitz gravity is re- 
stricted to a very tight window around its IR value. 

Finally, we should mention that although the present 
analysis provides the bounds for the running parameter, 
it does not enlighten the discussion about the possible 
conceptual problems and instabilities of Hofava-Lifshitz 
gravity, nor it can address the questions concerning the 
validity of its theoretical background, which is the sub- 
ject of interest of other studies. The present work just 
faces the problem from the phenomenological point of 
view, which is necessary but not sufficient, and thus 
its results can been taken into account only if Hofava- 
Lifshitz gravity passes successfully the aforementioned 
theoretical tests. 
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